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摘要 本 研究 针对 国产 高 剂量 率 近 距离 治疗 Ir-192 放射 源 ， 参 考 美 国医 学 物理 师 协会 AAPM) 在 TG43-U1 
中 推荐 的 剂量 学 参数 ， 采 用 蒙特 卡 罗 模 拟 软 件 建立 了 精细 结构 模型 ， 并 进行 了 模拟 计算 。 模 拟 计算 涵盖 了 剂 
量 率 常 数 A、 单 位 活 度 空气 比 释 动能 强度 S g [A 。”、 径 向 剂量 函数 g , r ”) 以 及 各 向 异性 函数 
F (r ,0 ). 结果 显示 ， 剂 量 率 常 数 计算 值 为 1.105 cGy. h' 。U1!， 与 文献 报道 的 差异 在 1. 2% 以 内 ; 
单位 活 度 空气 比 释 动能 率 为 9.788 X10 ”UBq ， 与 文献 报道 结果 的 差异 为 0. 23%。 研 究 还 获得 了 放射 源 中 垂 线 
0.5 至 20 厘米 范围 内 的 径 向 剂量 函数 ， 并 拟 合 得 到 了 经 验 公 式 。 整 体 而 言 ， 基 于 蒙特 卡 罗 软 件 建立 的 国产 
Ir-192 放射 源 模型 在 剂量 学 参数 上 与 文献 数据 具有 和 较 好 一 致 性 , 为 国产 Ir-192 源 的 临床 实践 应 用 提供 了 指导 。 
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Abstract [Background]: The Ir-192 high-dose-rate brachytherapy radiation source possesses the characteristic 
of high central dose rate with a rapid dose fall-off at the edges. This feature is highly advantageous in clinical 
applications, as it ensures that the tumor area receives sufficient radiation dose while significantly reducing the risk of 
damage to surrounding normal tissues and organs, thereby effectively protecting these healthy tissues. [Purpose]: 
To establish a detailed structural model of a domestically produced high-dose-rate Ir-192 brachytherapy source using 
Monte Carlo simulation software, based on the dosimetric parameters recommended by the American Association of 
Physicists in Medicine (AAPM) in TG43-UI, and to perform simulation calculations. [Methods]: Using the 
Monte Carlo software, a detailed structural model of the domestically produced high-dose-rate Ir-192 brachytherapy 
source was established. The simulated dosimetric parameters included the dose rate constant A, air kerma rate per 
unit activity S K i A „radial dose function g L (r  ),and anisotropy funcion F (r ,0 ). 
[Results]: The simulated dose rate constant was 1.105 cGy:h!:U-!, with a difference of less than 1.2% from the 
literature values. Air kerma rate per unit activity was 9.788 X 108 U- Bq'', with a difference of 0.2396 from the 


literature values. The radial dose function was obtained for distances from 0.5 to 20 cm from the source axis, and an 
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empirical formula was fitted. [Conclusions]: The domestically produced Ir-192 source model established using the 
Monte Carlo software shows good consistency with the literature-reported dosimetric parameters, indicating that this 
model can be used for clinical practice applications of domestically produced Ir-192 sources and has certain guiding 


significance. 
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Ir-192 近 距 治疗 源 是 一 种 高 剂量 率 Y 射线 放射 源 ， 其 具有 中 心 剂 量 率 高 且 边 缘 剂 量 迅 速 
下 降 的 特性 。 在 临床 治疗 过 程 中 , 这 种 剂量 分 布 特性 使 得 17-192. 放射 源 能 够 有 效 地 保护 肿瘤 
周围 的 正常 组 织 和 器 官 。 通 过 精确 控制 放射 源 的 位 置 和 照射 时 间 ， 能 够 实现 对 肿瘤 区 域 的 高 
效 照射 ， 同 时 将 对 周围 健康 组 织 的 损伤 降 至 最 低 。Ir-192 放射 源 因 
前 列 腺 癌 等 疾病 的 治疗 [1][2] 。 在 近 距 治疗 过 程 中 ， 剂 量 的 精准 计算 和 准确 分 布 直接 影响 患 
者 的 治疗 效果 。 为 了 确保 治疗 的 有 效 性 和 安全 性 , 必须 对 11-192. 放射 源 的 剂量 学 特性 进行 详 
细 的 定量 评估 。 美 国医 学 物理 学 家 协会 (AAPM) 在 其 发 布 的 TG-43[3] 和 TG-43U1[4] 报 告 
中 明确 建议 ， 在 临床 实施 前 ， 应 至 少 进行 一 套 完整 的 实验 测量 和 一 套 完整 的 蒙特 卡 罗 模 拟 ， 
以 获得 放射 源 的 剂量 学 参数 。 这 些 参数 对 于 治疗 计划 的 制定 和 实施 至 关 重 要 。 尽 管 国际 上 对 
Ir-192 放射 源 的 剂量 学 参数 已 有 大 量 研究 ， 但 国内 针对 原子 高 科 股份 有 限 公司 生产 的 HDR 
Ir-192 近 距 治疗 源 的 剂量 学 参数 的 模拟 研究 仍 较为 稀少 , 仅 有 吴 晓 维 等 人 研究 过 [5]。 为 了 使 
研究 数据 更 加 丰富 , 本 文采 用 蒙特 卡 罗 模 拟 方法 , 对 该 放射 源 的 剂量 学 参数 进行 了 系统 研究 。 
通过 使 用 蒙特 卡 罗 模 拟 软 件 , 建立 了 放射 源 的 精细 几何 模型 , 分 别 计算 光子 在 空气 中 的 空气 
比 释 动能 强度 和 水 介质 中 的 吸收 剂量 。 随 后 ， 将 模拟 结果 与 现 有 文献 数据 [10-16] 进 行 对 比 ， 
以 验证 本 研究 的 准确 性 和 可 靠 性 ,这 一 工作 为 进一步 提高 近 距 离 放射 治疗 的 剂量 计算 精度 提 
供 了 重要 的 理论 依据 和 数据 支持 。 
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1 材料 与 方法 
1. 原子 高 科 HDR Ir-192 源 结 构 

图 1 展示 了 由 原子 高 科 股份 有 限 公 司 生产 的 HDR Ir-192 放射 源 的 示意 图 。 根 据 制 造 商 
提供 的 数据 ， 该 放射 源 由 一 个 直径 为 0.6 毫米 、 长 度 为 3.5 毫米 的 圆柱 体 组 成 ， 采 用 纯 久 金 
属 制 成 ， 其 密度 为 22.56 g/cm? . Ir-192 是 一 种 半衰期 约 73.8 天 的 放射 性 同位 素 ， 平 均 每 次 
衰变 会 释放 出 一 个 电子 和 2.36 个 光子 ， 其 放射 性 物质 在 圆柱 体内 均匀 分 布 。 该 放射 源 被 封 
装 在 一 个 长 6.5 毫米 的 06Cr19Nil0 (304 不 锈 钢 ) 圆柱 体 中 ， 圆 柱 的 内 径 为 0.8 毫米 ， 外 径 
为 1.1 毫米 。 放 射 源 一 端 连接 一 根 长 度 为 170 厘米 的 06Cr19NilO (304 不 锈 钢 ) W, IAk 
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只 展示 了 6cm; 另 一 端 为 半径 0.55 毫米 的 半圆 球状 结构 ， 在 放射 源 与 封装 包 壳 之 间 的 空间 
中 充满 空气 。 封 装 的 06Crl19Nil10〈304 不 锈 钢 ) 具有 优良 的 机 械 强 度 和 耐 腐蚀 性 ， 能 够 有 效 
地 保护 内 部 的 放射 性 物质 ， 并 确保 其 在 使 用 过 程 中 不 发 生 泄漏 。 表 1 详细 列 出 了 在 计算 模 
型 中 使 用 的 各 种 材料 的 组 成 和 密度 [6][7]， 为 后 续 的 蒙特 卡 罗 模 拟 计算 提供 了 基础 数据 。 这 
些 参数 对 于 准确 模拟 放射 源 的 剂量 分 布 至 关 重 要 ， 确 保 了 计算 结果 的 可 靠 性 和 有 效 性 。 
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图 1 原子 高 科 HDR Ir-192 源 结 构 示意 图 
Fig.1 Schematic diagram of Atomic Hi-Tech Co., Ltd. HDR Ir-192 source structure 


表 1 物质 的 材料 组 成 及 百分比 


Table1 Composition and Percentage of Materials 


材料 密度 元 素 及 成 分 

Material Density / g:cm ? Elements and ingredients 

久 Iridium 22.56 Ir(1) 

06Cr19Ni10 793 C(0.08). P(0.045). Si(1.0). S(0.03). 


Ci(19), Mn(2.0), Fe(68.345). Ni(9.5) 
H(0.0124). N(23.1781). O(75.5267). 
Ar(1.2827) 

液态 水 Water 0.998207 H(11.1894). O(88.8106) 


e (TM 1 
ZA CFH) Air (dry) 1205 x 103 


2. 剂量 学 参数 计算 公式 


美国 医学 物理 师 协会 (American Association of Physicists in Medicine,AAPM) 建议 , 在 将 


每 种 新 近 距 离 放射 治疗 源 用 于 治疗 计划 之 前 ， 必 须 依 据 TG43-U1 报告 计算 其 剂量 学 参数 [4]。 


根据 AAPM 的 TG-43U1 报告 ， 这 些 剂量 学 参数 包括 : 剂量 率 常数 A， 单 位 活 度 空气 比 释 动 


能 强度 S K /A ”， 径 向 剂量 函数 g , (rr ”)， 各 向 异性 函数 F (r ,0 ). 遵循 


TG43 报告 建议 的 剂量 计算 规范 ， 在 介质 中 某 点 (r ,0 ) 的 吸收 剂量 率 D (r ,0 ) 


在 极 坐标 下 表 式 为 : 


D (r 8 )=S pA :一 上 


gJ L (r F (r ,0 ) (1) 


其 中 ，7 是 距 活性 源 中 心 的 径 向 距离 ，9 ”是 相对 于 源 的 纵 轴 的 极 角 ,， G , (r 


,0 ) 是 几何 因子 ，S k 是 放射 源 的 空气 比 释 动能 强度 ， 单 位 为 U 或 者 cGy:cm2.h1， 其 


计算 公式 为 : 
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o (d E (2) 


RP, K o (d ) 为 在 真空 下 自由 空间 中 ， 放 射 源 横 轴 上 d om 处 的 空气 比 释 动能 率 ， 


单位 为 cGy:h"!; 4 为 探测 点 到 放射 源 横 轴 中 心 的 距离 ， 为 100cm。 


"D r ,6 ) 和 5S 


k 已 知 时 , 可 以 通过 公式 计算 A ig 


Lo DEF 


(r ,0 ”)。 剂 量 率 常数 A 是 指 在 参考 点 位 置 (rY ” 0,0 0) 处 的 水 吸收 剂量 率 与 空气 比 
释 动能 强度 的 比值 ， 其 计算 公式 如 下 : 


jh eT oo) 3) 
S K 


xu, D (r 909r 0=lem，9 ”0=90" 位 置 处 的 吸收 剂量 率 ，A ”的 单位 为 


cGy-h U^; 


为 径 向 剂量 函数 ， 为 各 向 异性 函数 ， 其 计算 公式 如 下 : 


D (r ,0 0) G L (r 0,0 0) 
r = 一 d (4) 
gr ) D r 00 0) CLT ,0 0) 


Der ,0g )G Ylr ,0 o9 
F (r ,0 二 一 一 一 一 一 :一 一 一 一 一 一 (5) 
( ) Dr ,9 0 Gr (lr ,0) 


AF, G , 上 0) 为 线 源 几 何 因 子 ， 其 数值 取决 于 放射 源 的 几何 尺寸 ， 表 示 源 自身 对 其 周 


围 剂量 率 的 影响 。 该 几何 因子 的 计算 公式 为 : 


人 工 (r ,0 Js T "ET 0 +O 
G 1 C ,6 )=(r ?-L "METEO 0 =0 


式 中 ， 世 为 放射 源 有 效 活 性 长 度 ， 有 ”为 任意 一 点 对 放射 源 两 端 夹 
3. 蒙特 卡 罗 模 拟 


本 研究 采用 蒙特 卡 罗 模 拟 软件 ， 为 了 更 方便 、 准 


zu 


描述 放射 源 的 几何 结构 和 模拟 计算 ， 
在 极 坐 标 系 下 ,建立 了 用 于 计算 Ir-192 放射 源 剂量 学 参数 的 模型 ， 如 图 2 所 示 。 光 子 截面 数 
据 库 选用 McpLib04, 以 确保 模拟 的 准确 1 


ms 


生 。 通 过 源 卡 对 放射 源 进行 设 定 ， 其 部 分 程序 代码 为 : 


SDEF POS--0.195 0 0 ERG=D1 PAR-2 AXS-1 0 0 RAD-D2 EXT=D3。 最 后 计算 了 穿 过 放射 源 包 
壳 的 空气 比 释 动能 强度 以 及 水 介质 中 的 吸收 剂量 。 
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06Cr19Ni10 Ar B = 


图 2 剂量 参数 的 参考 坐标 系 


Fig.2 Reference coordinate system for dose parameters 


水 吸收 剂量 的 计算 : 将 放射 源 模型 放置 于 一 个 的 长 50 cm, #4 25 cm 的 圆柱 形 水 模 体 


中 心 处 ， 并 假设 水 组 织 内 带电 粒子 处 于 平衡 状态 。 使 用 点 探测 器 的 F4 计数 卡 和 NIST 官网 


XCOM 数据 库 中 的 水 的 质量 
置 的 水 的 吸收 剂量 。 在 计算 径 向 剂 


影响 ， 我 们 采用 在 不 同 距离 位 置 划分 不 同 的 体 元 的 策略 : 在 距离 中 心 + X Sem 时 ， 采 用 半 


an 
CC 


量 吸收 系数 ， 或 者 使 用 探测 器 的 F6 能 量 沉积 卡 ， 记 录 不 同位 


fg 7 (7 ”) 时 , 为 尽 可 能 减少 体 元 大 小 对 结果 的 


径 为 0.035cm 的 球形 体 元 ; 在 5cm <r < 9 em 范围 内 ， 采 用 半径 为 0.055cm 的 球形 体 元 ; 


在 9cm<r < 15cm 范围 内 ， 采 用 半径 为 0.06cm 的 球形 体 元 ; 在 15 cm<r < 20 em 范围 


内 ,采用 


EIEN 0.075cm 的 球形 体 元 [12] .为 了 减少 体 元 大 小 对 各 向 异性 函数 F (r ,0 ) 


模拟 过 程 中 角度 的 影响 , 本 研究 在 不 同 径 向 距离 x 和 0” 一 180” 和 角度 范围 内 , 采用 F5 卡 (点 
探测 器 ) 进行 模拟 。 模 拟 过 程 中 设置 光子 和 电子 的 输 运 模式 〈mode: p,e) ， 考 虑 初级 光子 和 
次 级 电子 的 贡献 。 粒 子 数 设 定 为 1X108 ， 光 子 和 电子 的 截止 能 量 分 别 设 定 为 10 keV 和 2 


MeV[13]。 


2.1 剂量 率 常数 和 单位 活 度 空气 比 释 动 能 强度 : 利用 蒙特 卡 罗 软 件 分 别 计算 了 , 在 真空 环 


境 下 距 源 径 向 距离 为 100cm 处 的 空气 比 释 动能 强度 , 以 及 在 水 的 环境 下 距 源 径 向 lcm 处 的 吸 
收 剂量 率 , 经 过 公式 计算 及 单位 转化 ,本 研究 得 出 剂量 率 常 数 为 1.105cGy:h-1.U!, 与 文献 [14] 


中 的 1.113cGy.h1.U-1 偏差 在 0.72%; 与 文献 [10,15] 中 的 1.115 cGy-ht-U7 FI 1.117 eGy-hi-U" 


1 


扁 差 在 0.90% 和 1.09%; Xi 


Zu 
AUTE 


导出 单位 活 度 空 气 比 释 动能 强度 为 9.788 X 108 U-Bq'' 和 文献 


[10,15] 中 的 9.781X103U.Bq1 和 9.790 X 105 U-Bq'! 偏差 在 0.069% 和 0.023%. 
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2.2 径 向 剂量 函数 : 基于 第 二 节 建 立 的 模拟 环境 , 计算 得 到 了 放射 源 在 不 同 径 向 位 置 的 数 


据 。 为 了 提高 模拟 结果 的 准确 性 和 精细 度 ， 我 们 在 径 向 方向 上 增加 了 探测 器 的 数量 。 将 模拟 
所 得 数据 代入 公式 〈4) ， 计 算得 出 随 着 放射 源 径 向 距离 增加 而 剂量 率 降低 的 径 向 剂量 函数 。 
随后 ， 将 本 研究 所 得 的 径 向 剂量 函数 与 相关 文献 [13,16] 中 的 研究 数据 进行 了 对 比 ， 整 理 结果 
如 图 3 所 示 。 
在 径 向 距离 20cm 以 内 ， 将 本 研究 结果 与 相关 文献 [13] 相 比 ， 相 对 误差 范围 为 0.26% 一 
3.76%; 与 文献 [16] 相 比 ， 相 对 误差 范围 为 0.27% 一 3.79%， 显 示 出 较 好 的 一 致 性 。 在 径 向 距 
离 r=0.05cm 至 5cm 范围 内 ， 径 向 剂量 函数 呈 上 升 趋势 ， 与 文献 [13] 的 报道 基本 一 致 : 在 
径 向 距离 3 cm 以 内 ， 径 向 剂量 函数 出 现 小 幅 波动 ， 表 现 为 先 下 降 升 后 上 升 的 趋势 。 当 径 向 距 
离 r«6cm 时 ， 所 有 研究 结果 显示 镀 源 的 径 向 剂量 函数 近似 保持 不 变 ， 而 当 径 向 半径 超过 6 
cm 时 ， 径 向 剂量 函数 逐渐 下 降 。 当 径 向 半径 超过 15 cm 后 ， 本 研究 的 径 向 剂量 函数 与 文献 
[13,16] 的 结果 差异 逐渐 增 大 。 

本 研究 利用 模拟 获得 的 数据 ， 通 过 Origin 绘图 工具 对 径 向 剂量 函数 关于 距离 x 进行 五 阶 


IT 


多 项 式 拟 合 [17]， 得 出 了 以 下 经 验 公 式 : g , C )-(a ota 1xr +a ax 


r ?+a 3xr 3+a axr 4+a sxr 5， 在 径 向 距离 r<20cm 的 范围 内 ， 


拟 合 系数 分 别 为 0.98206. 0.00877. - $.31377X104、- 1.72573 X 104、1.14432X105 和 - 


2.37093 X 107。 拟 合 图 形 如 图 4 所 示 。 


=— This study 
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^  GammaMed 


Radial dose function g,(r) 
x 
"» 


T T T T T T T T 
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图 3 径 向 剂量 函数 对 比 图 


Fig.3 Radial dose function comparison Figure 
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图 4 径 向 剂量 函数 拟 合 曲线 
Fig.4 Fitting curve of radial dose function for HDR 192-Ir source 


3. 各 向 异性 函数 : 在 本 研究 中 ,根据 探测 器 


E 放 射 源 不 同 昌 


E 离 (0.5cm. lcm, 2cm, 3cm, 


4cm, 5cm, 6cm, 7cm, 8cm, 9cm, 10cm, 15cm, 20cm) 、 不 同 角度 (0° ~180° ) 得 到 


模拟 数据 ， 带 入 到 公式 〈5) 中 ， 得 到 二 维 各 向 异性 函数 ， 并 将 其 整理 于 表 2 中 ， 将 表 2 中 


lcm, 3cm, 5cm, 7em 不 同 径 向 距离 下 得 到 的 各 向 异性 函 与 文献 [13,16] 对 比 ， 绘 制 成 如 图 5 


所 示 图 像 。 在 将 本 研究 与 文献 [13,16] 中 的 各 向 异性 函数 对 比 ， 通 过 对 比分 析 可 以 发 现 ， 当 


r=10cm, 0 =180° 时， 各 向 异性 函数 函数 偏差 明显 1 


源 模 型 的 源 辫 结构 不 同 所 导致 的 ， 在 180” 以 内 


误差 为 3.35%， 所 得 到 的 模拟 结果 具有 较 好 的 一 致 性 。 


WX, dx 


是 由 于 模拟 过 程 中 所 建立 放射 


， 本 研究 成 果 


F 


与 文献 [13,16] 对 比 ， 最 大 相对 


X2 HDR 192-Ir 各 向 异性 函数 表 


Table2 Table of anisotropy function values for HDR 192-Ir source 


0 0. 5cm lem 2cm 3cm 4cm 5cm 6cm Tcm 8cm 9cm 10cm 15cm 20cm 

j 0.6533 0.6190 0.6131 0.6364 0.6559 0.6808 0.7017 0.7211 0.7311 0.7519 0.7604 0.8018 0.8295 
10* 0.7194 0.7085 0.7251 0.7427 0.7572 0.7734 0.7877 | 0.8018 0.8066 0.8180 0.8221 0.8588 0.8798 
20° 0.8440 0.8332 0.8395 0.8454 0.8544 0.8659 0.8740 0.8761 0.8807 0.8853 0.8875 0.9053 0.9221 
30* 0.9084 0.9004 0.9040  À 0.9088 0.9093 0.9138 0.9233 0.9250 0.9251 0.9323 0.9319 0.9442 0.9594 
40° 0.9462 0.9414 0.9447 0.9446 0.9449 0.9481 0.9521 0.9569 0.9499 0.9530 0.9572 0.9636 0.9857 
50* 0.9689 0.9678 0.9684 0.9688 0.9712 0.9702 0.9778 0.9758 0.9711 0.9817 0.9802 0.9767 1.0081 
60° 0.9832 0.9833 0.9829 0.9849 0.9822 0.9858 0.9883 0.9876 0.9849 0.9856 0.9894 0.9828 0.9997 
70° 0.9923 0.9938 0.9954 0.9942 0.9911 0.9938 0.9989 1.0031 0.9953 0.9941 1.0040 0.9907 1.0115 
80° 0.9982 0.9984 0.9988 1.0008 0.9982 1.0004 1.0024 1.0015 1.0028 1.0004 1.0006 1.0337 0.9941 
90* 1.0000 1.0000 1.0000  Á 1.0000  Á 1.0000  Á 1.0000 1.0000  Á 1.0000  Á 1.0000 1.0000  Á 1.0000  Á 1.0000 1.0000 
100* 0.9989 0.9968 1.0059  Á 0.9980 0.9975 0.9979 0.9991 1.0006 0.9964 1.0105 0.9991 0.9917 1.0056 
110* 0.9922 0.9932 0.9970 0.9938 0.9886 0.9978 0.9966 0.9991 1.0125 0.9962 0.9960 0.9946 0.9934 
120* 0.9839 0.9839 0.9832 0.9858 0.9809 0.9877 0.9893 0.9953 0.9843 0.9942 0.9842 0.9839 0.9921 
130* 0.9689 0.9684 0.9684 0.9697 0.9666 0.9705 0.9742 0.9802 0.9766 0.9813 0.9719 0.9743 0.9922 


Eh 技 术 20XX, XX: XXXXXX 


140* 0.9458 | 0.9417 0.9489 0.9462 0.9487 0.9480 0.9537 0.9539  Á 0.9547 0.9598 0.9563  Á 0.9573 0.9806 
150* 0.9018 | 0.9032 0.9071 0.9079 0.9172 0.9188 0.9250 0.9225 0.9277 0.9334 0.9356 0.9575 
160* 0.8330 0.8384  Á 0.8484  À 0.8527 0.8605 0.8707 0.8838  Á 0.8787 0.8863 0.8870 0.9085 0.9190 
170* 0.7125 0.7318 0.7443 0.7596 0.7771 0.7872 0.7973 | 0.8084  À 0.8151 0.8460 0.8614 
180* 0.4947 0.5243 0.5541 0.6572 0.7193 
w— This study 
*  GammaMed Plus 
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图 4 本 研究 和 文献 报道 的 距 源 中 心 不 同 距离 的 各 向 异性 函数 对 比 图 


Fig.4 Comparison of the anisotropy function at different distances from the source center between this study and the literature 


reports 


3 讨论 


LT 


本 研究 利用 蒙特 卡 罗 模 拟 软件 ， 在 极 坐 标 系 下 建立 了 用 于 计算 1-192 放射 源 剂量 学 参数 
的 模型 。 根 据 TG-43U1 研究 报告 中 的 剂量 计算 公式 ， 模 拟 计 算 了 原子 高 科 股份 有 限 公司 
HDRIr-192 源 的 剂量 学 参数 。 首 先 ， 模 拟 计算 了 单位 活 度 空气 比 释 动能 强度 ， 结 果 显 示 本 帮 
究 得 出 的 剂量 率 常 数 为 1.105cGyh-IU-: 与 文献 [14] 中 的 1.113cGyh-1U-! 相 比 ， 偏 差 为 0.72%; 与 


文献 [10,15] 中 的 1.115 cGyh U^ 和 1.117 cGyh-1U-! 相 比 , 偏差 分 别 为 0.90% 和 1.09%， 考 虑 放 
射 源 结构 与 文献 中 放射 源 结构 上 的 偏差 ， 结 果 偏 差 范 围 可 接受 。 其 次 ， 本 文 计算 了 径 向 剂量 
函数 , 并 使 用 origin 绘图 工具 对 径 向 剂量 函数 关于 距离 ,进行 了 五 阶 多 项 式 拟 合 , 与 文献 [13,16] 
的 结果 进行 比较 ， 发 现 径 向 剂量 函数 在 整体 趋势 上 基本 一 致 ， 但 随 着 半径 的 增 大 ， 模 拟 结果 
的 偏差 也 有 所 增加 ， 这 可 能 是 由 于 放射 源 几何 结构 的 不 同 所 致 。 最 后 ， 本 文 研究 了 各 向 异性 
函数 ， 通 过 比较 发 现 ， 本 文 所 得 数据 与 已 发 表 文 献 的 结果 在 误差 范围 内 保持 一 致 ， 验 证 了 模 
拟 方法 的 可 靠 性 和 结果 的 准确 性 。 本 研究 为 国产 Ir-192 放射 源 提供 了 一 种 可 靠 的 剂量 率 常 数 
计算 方法 , 给 国产 放射 源 剂量 学 进展 提供 了 数据 文 撑 , 以 提高 近 距 离 治 疗 临 床 治疗 的 准确 性 。 
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